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Abstract 



We revisit a class of supersymmetric SO(10) models with t-b-r Yukawa coupling 
unification condition, with emphasis on the prediction of the Higgs mass. We discuss 
!>> qualitative features in this model that lead to a Higgs mass prediction close to 125 

GeV. We show this with two distinct computing packages, Isajet and SuSpect, and 
also show that they yield similar global features in the parameter space of this model. 
We find that t-b-r Yukawa coupling unification prefers values of the CP-odd Higgs 
mass rriA to be around 600 GeV, with all colored sparticle masses above 3 TeV. We 
also briefly discuss prospects for testing this scenario with the ongoing and planned 
r^i direct dark matter detection experiments. In this class of models with t-b-r Yukawa 

unification, the neutralino dark matter particle is heavy (rn^o > 400 GeV), which 
coannihilates with a stau to yield the correct relic abundance. 
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1 Introduction 



Despite no signals from direct Supersymmetry (SUSY) searches at the LHC, a 125 
GeV Higgs boson is a strong indication that a SUSY signal might be imminent. SUSY 
still remains at the forefront of beyond the standard model physics scenarios due to 
several reasons. In addition to solving the gauge hierarchy problem and providing 
a dark matter candidate, it also leads to unification of the gauge couplings. In the 
SUSY case, t-b-r Yukawa unification (YU) [1] can also be accommodated in contrast 
to its non-SUSY version. Both these observations hint at an underlying grand unified 
structure like SO (10) which may be supersymmetric. The implications of t-b-r YU 
has been extensively explored over the years [2, 3]. 

The discovery of the Higgs boson was announced by the ATLAS and the CMS 
collaborations using the combined 7 TeV and 8 TeV data, with ATLAS observing 
a 5.0(j signal at 126.0 ± 0.4(stat) ± 0.4(syst) GeV [4], and CMS a 5.0a signal at 
125.3 ± 0.4(stat) ± 0.5(syst) GeV [5]. This observation was also confirmed at the 
de Moriond 2013 conference. Spurred by these exciting results, we revisit a previous 
analysis done in ref. [6] with special emphasis on the prediction of the Higgs mass in 
a class of SO(10) models with t-b-r YU condition. It was shown in [6] that with non- 
universal soft supersymmety breaking (SSB) gaugino masses, which can be derived in 
the framework of SO(10) GUT, and with universal SSB scalar Higgs doublet masses 
at M GUT (Mh u = Mff ), the mass of CP-even SM-like Higgs boson can be predicted 
with t-b-r YU case. In this paper we shed light on the qualitative features that lead to 
consistency of t-b-r YU with a 125 GeV Higgs. Furthermore, we employ two different 
computing packages, namely Isajet 7.84 and SuSpect 2.41, and show that they agree 
very well with the 125 GeV Higgs mass prediction for t-b-r YU. Good agreement is 
found between these two programs over the entire parameter space, thereby rendering 
our conclusions more robust. 

Another motivation to revisit the analysis presented in ref. [6] is the recent dis- 
covery of B s — > decay by the LHCb collaboration [7]. The branching frac- 
tion BF(B S -> fi + /J,-) = 3.2±H x 10~ 9 is in accord with the SM prediction of 
(3.2 ±0.2) x 10~ 9 [8]. In SUSY models, this flavor- changing decay receives contribu- 
tions from the exchange of the pseudoscalar Higgs boson A [9] , which is proportional 
to (tan f3) 6 /m\. Since tan (5 « 47 was predicted in ref. [6], it is interesting to see how 
the parameter space is impacted with the B s — > discovery. 

We also intend to highlight another interesting feature of this YU model, namely 
that in addition to the prediction of a CP even Higgs boson mass of around 125 
GeV , the model also prefers the CP odd Higgs boson mass of around 600 GeV. 
This prediction can hopefully be tested at the LHC in the near future [10]. The 
colored sparticle masses, consistent with good (10% or better) t-b-r YU, lie well 
above the current mass limits from the LHC, i.e., m g > 1.4 TeV (for m- g ~ m,) and 
nig > 0.9 TeV (for rrig < m$)[ll, 12]. 
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The outline for the rest of the paper is as follows. In section 2 we briefly de- 
scribe an SO(10) GUT model which can yield realistic fermion masses and mixings 
compatible with t-b-r YU. In Section 3 we present the parameter space that we scan 
over, and describe how the MSSM gaugino mass relations can be obtained at Mgut- 
In Section 4 we summarize the scanning procedure and the experimental constraints 
applied in our analysis. In section 5 we discuss how SUSY threshold corrections, 
which are necessary to obtain t-b-r YU and radiative electroweak symmetry breaking 
(REWSB), determine the CP-even Higgs boson mass. In section 6 we compare the 
results from the two packages, Isajet and SuSpect. In section 7 we briefly discuss the 
little hierarchy problem in the presence of t-b-r YU, and our conclusions are presented 
in Section 8. 



2 SO (10) GUT with t-b-r Yukawa unification 



One of the main motivations of SO(10) GUT, in addition to gauge coupling unifi- 
cation, is matter unification. The spinor representation of SO(10) unifies all matter 
fermions of a given family in a single multiplet (16»), which also contains the right 
handed neutrino (i/r) that helps to generate light neutrino masses via the see-saw 
mechanism [13] . The right handed neutrino can also naturally account for the baryon 
asymmetry of the Universe via leptogenesis [14]. Another virtue of SO(10) is that, 
in principle, the two MSSM Higgs doublets can be accommodated in a single ten 
dimensional representation, which then yields the following Yukawa couplings 



Here i,j = 1,2,3 stand for family indices and the SO(10) indices have been omitted 
for simplicity. For the third generation quarks and leptons, the interaction in Eq.(l) 
yields the following Yukawa coupling unification condition at the GUT scale [1] 



It is interesting to note that in the gravity mediation SUSY breaking scenario [15], 
t-b-r YU condition leads to LHC testable sparticle spectrum [2, 3] and even predicts 
a 125 GeV light CP-even Higgs boson mass [6]. On the other hand, it is well known 
that the interaction in Eq.(l) leads to a naive SO(10) relation: N = U oc D — L, 
where U, D, N and L denote the Dirac mass matrices for up and down quarks, 
neutrinos and charge leptons respectively. U oc D would imply vanishing quark flavor 
(Cabibbo-Kobayashi-Maskawa (CKM)) mixing [16], and m[!/m t = m®/m® which 
significantly contradicts with experimental observations. The superscript zero refers 
to the parameters evaluated at the Mgut- D = L T , which is a naive SU(5) relation, 
would imply m" = m° and m° d = m° e and this relation also strongly disagrees with 



Yij 16; 16; 10 H . 



(1) 



Y t = Y b = Y T = Y VT . 



(2) 
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the measurements. So, it is obvious the interaction in Eq.(l) should be modified in 
order to accommodate the observed pattern of quarks and mixing. 

There are two main approaches to avoid these problems and obtain realistic 
fermion masses and mixings. One way is to extend the Higgs sector and assume 
that the SM Higgs doublet fields are superposition of fields from the different SO(10) 
representations [17]. Another way is to introduce additional vector-like matter mul- 
tiplets at the GUT scale [18] which mix in a nontrivial way with fermions in the 16 
dimensional representation. It is equivalent to introducing non-renormalizable cou- 
plings which involve a non-singlet SO(10) field that develops a VEV [19]. Both these 
cases, however, allow t-b-r YU to be completely destroyed or partially destroyed to 
yield b-r YU. Therefore, in order to maintain t-b-r YU the following two conditions 
should be satisfied: 

• Third generation charged fermions obtain their masses only from Eq.(l). 

• The MSSM Higgs doublets H u and Hd reside solely in 10#. 

There exist models where t-b-r YU is maintained to a very good approximation with 
realistic fermion masses and mixings [20]. 

We consider the case in which the MSSM Higgs fields are contained in the 10 h 
representation and assume non-renormalizable couplings in the superpotential 

Y„ 16. 16, 10 4- r v 16, 16, (§)" 10 4- /, ^ ^ T& " . (3) 

Here M, a scale associated with the effective non-renormalizable interaction, could 
plausibly lie somewhere between the unification scale and the reduced Planck mass 
(Mp = 2.4 x 10 18 GeV). In Eq.(3) n is an integer and 45 is an adjoint representation of 
SO (10). We can have several 45-dimensional fields with VEVs in different directions 
in the space spanned by the 45 generator of SO(10), as long as the SM gauge group is 
unbroken. Ref. [19] shows that in order to have the correct naive SO(10) and SU(5) 
relationships for the first two families we need to have at least two 45 dimensional 
fields with VEVs in different directions in the SO(10) space. If one of the 45-plet 
develops a VEV along the B — L direction it can naturally lead to the so-called Georgi- 
Jarlskog relation [21], m° s = m°/3 and mJJ = ?>m° e . This will give correct masses after 
evolution of the RGEs to low scale. The second 45-plet pointing towards I 3 r, the 
third component of the right handed isospin group, will break U oc D and will allow 
for the correct CKM mixing matrix. 

The pair of spinors 16h + 16// break the rank of the group from five to four 
and provides a mass for the right handed neutrinos. The third term in Eq.(3) is 
responsible for generating the right handed neutrino masses. The adjoint 45-plet 
completes the breaking of SO(10) to the SM gauge group. In order to keep the 
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MSSM Higgs field in 10#, the term hi 16# 16# 10 allowed in principle, should not 
be in the superpotential. We can just assume the hi coupling to be zero, or it can 
be forbidden by introducing an additional symmetry. It was shown in [19] that the 
superpotential in Eq.(3) perfectly describes the observed fermion masses and mixings. 

Having the MSSM Higgs fields in the 10# representation, the doublet-triplet split- 
ting problem can be solved in SO(10) using the missing VEV mechanism [22]. A 
variety of realistic models based on this mechanism have been constructed [23]. 



3 Fundamental Parameter Space 

It has been pointed out [24] that non-universal MSSM gaugino masses at M GU t can 
arise from non-singlet F-terms, compatible with the underlying GUT symmetry such 
as SU(5) and SO(10). The SSB gaugino masses in supergravity [15] can arise, say, 
from the following dimension five operator: 

rpab 

\ a \ b + c.c. (4) 



2Mi 

Here A a is the two-component gaugino field, F ab denotes the F-component of the 
field which breaks SUSY, the indices a, h run over the adjoint representation of the 
gauge group. The resulting gaugino mass matrix is (F ab )/Mp, where the supersym- 
metry breaking parameter (F ab ) transforms as a singlet under the MSSM gauge group 
SU(3) C x SU{2)l x U(1)y- The F ab fields belong to an irreducible representation in 
the symmetric part of the direct product of the adjoint representation of the unified 
group. This is a supersymmetric generalization of operators considered a long time 
ago [25]. 

In SO(10), for example, 

(45 x 45) 5 = l + 54 + 210 + 770. (5) 

If F transforms as a 54 or 210 dimensional representation of SO(10) [24], one obtains 
the following relation among the MSSM gaugino masses at Mqut : 

M 3 : M 2 : M x = 2 : -3 : -1, (6) 

where Mi,M 2 ,M 3 denote the gaugino masses of U(l), SU(2)l and SU(3) C respec- 
tively. The low energy implications of this relation have recently been investigated 
in [26] without imposing YU. In this paper we consider the case with \x > and non- 
universal gaugino masses defined in Eq.(6). In order to obtain the correct sign for the 
desired contribution to (g — 2) M , we set Mi > 0, M 2 > and M 3 < 0. Somewhat to 
our surprise, we find that this class of t-b-r YU models make a rather sharp prediction 
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for the lightest SM-like Higgs boson mass [6] . In addition, lower mass bounds on the 
masses of the squarks and gluino are obtained. 

Notice that in general, if F ab transforms non trivially under SO(10), the SSB 
terms such as the trilinear couplings and scalar mass terms are not necessarily uni- 
versal at M GUT . However, we can assume, consistent with SO(10) gauge symmetry, 
that the coefficients associated with terms that violate the SO(10)-invariant form are 
suitably small, except for the gaugino term in Eq.(6). We also assume that D-term 
contributions to the SSB terms are much smaller compared with contributions from 
fields with non-zero auxiliary F-terms. 

Employing the boundary condition from Eq.(6), one can define the MSSM gaugino 
masses at Mqut i n terms of the mass parameter M1/2 : 

M x = Mi/a 
M 2 = 3Mi /2 

M 3 = -2M 1/2 . (7) 

Note that M 2 and M 3 have opposite signs which, as we will show, is important in 
implementing Yukawa coupling unification to a high accuracy In order to quantify 
Yukawa coupling unification, we define the quantity R t br as, 

max(yt,y b ,y T ) , Q , 
mm(yt, y b , y T ) 

We have performed random scans for the following parameter range: 

< m 16 < 10 TeV 
< m 10 < 10 TeV 
< Mi /2 < 5 TeV 
35 < tan/3 < 55 
-3 < A /m 16 < 3 (9) 

Here rriiQ is the universal SSB mass for MSSM sfermions, mio is the universal SSB 
mass term for up and down MSSM Higgs masses, M1/2 is the gaugino mass parameter, 
tan (3 is the ratio of the vacuum expectation values (VEVs) of the two MSSM Higgs 
doublets, Aq is the universal SSB trilinear scalar interaction (with corresponding 
Yukawa coupling factored out). We use the central value m t = 173.1 GeV and lcr 
deviation (m t = 174.2 GeV) for top quark in our analysis [27]. We choose a +lcr 
deviation in m t since it leads to an increase in the Higgs mass and improves the 
prediction of the Higgs mass in our analysis. Our results however are not too sensitive 
to one or two sigma variation in the value of m t [28]. We use mb(mz) = 2.83 GeV 
which is hard-coded into Isajet. 
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4 Constraints and Scanning Procedure 



We employ Isajet 7.84 [29] and SuSpect 2.41 [30] interfaced with Micromegas 2.4 [31] 
to perform random scans over the fundamental parameter space. Isajet and SuSpect 
employ full two loop RGEs for the SSB parameters between M z and M GUT . The 
approach employed by both is similar, but there are some important differences that 
have been previously studied [32, 33]. SuSpect assumes that the full set of MSSM 
RGEs are valid between Mz and Mqut and uses what is referred to as the 'common 
scale approach' in [33]. In this approach the DR parameters are extracted at a 
common scale M EWS b = yj m i 1 m i 2 - Isajet, on the other hand, uses the 'step-beta 
function approach' which means that it employs one-loop step beta functions for 
gauge and Yukawa couplings. These two approaches yield very similar results for 
most of the SUSY parameter space. There are some regions where the discrepancies 
get magnified [32] and which we observe in our results as well. These differences will 
be discussed in section 6. 

An approximate error of around 2 GeV in the estimate of the Higgs mass in Isajet 
and SuSpect largely arise from theoretical uncertainties in the calculation of the min- 
imum of the scalar potential, and to a lesser extent from experimental uncertainties 
in the values for m t and a s . 

An important constraint comes from limits on the cosmological abundance of sta- 
ble charged particles [35] . This excludes regions in the parameter space where charged 
SUSY particles become the lightest supersymmetric particle (LSP). We accept only 
those solutions for which one of the neutralinos is the LSP and saturates the WMAP 
bound on the relic dark matter abundance. 

Micromegas is interfaced with Isajet and SuSpect to calculate the relic density and 
branching ratios BR(B S — > fi + n~) and BR(b — > 57). With these codes we implement 
the following random scanning procedure: A uniform and logarithmic distribution of 
random points is first generated in the parameter space given in Eq. (9). The func- 
tion RNORMX [36] is then employed to generate a gaussian distribution around each 
point in the parameter space. The data points collected all satisfy the requirement 
of radiative electroweak symmetry breaking (REWSB), with the neutralino in each 
case being the LSP. After collecting the data, we impose the mass bounds on all the 
particles [35] and use the IsaTools package [37] to implement the various phenomeno- 
logical constraints. We successively apply the following experimental constraints on 
the data that we acquire from SuSpect and Isajet: 

< 6.2 x 1(T 9 (2a) [38] 

< 3.87 x 1(T 4 (2a) [39] 

< 2.41 (3a) [39] 

< 55.6 x 1(T 10 [40] 



0.8 x 10~ 9 < BR(B S ->- 

2.99 x 10- 4 < BR(b si) 

Q 15 < BR(B u ^>-ti/ t )mssm 
— BR(B u ^tru T )sM 

< A(g - 2)J2 
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mis ( TeV) mis ( TeV) 

(a) Isajet 7.84 (b) SuSpect 2.41 

Figure 1: Plot in the M 3 — m 16 planes. The panel on the left shows data points 
collected with Isajet 7.84 whereas the panel on the right shows data from SuSpect 
2.41. The light blue points are consistent with REWSB and neutralino LSP. For both 
the panels, black points are subset of the light blue points and satisfy R t b T < 1-03 in 
panel 1(a) and R t b T < 1-05 in panel 1(b). The unit line is to guide the eye. 
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-1 0.5 0.5 1 

A tJ V m t L m iR 

Figure 2: Plot in the J m lt m tR ~~ A/ yJ m t L m t R planes. The data points shown are 
collected using Isajet. Color coding is the same as in Figure 1. 
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(c) Isajet 7.84 + SuSpect 2.41 



Figure 3: Plot in the Rtb T — n^h plane. Panel 3(a) shows results obtained from Isajet, 
while 3(b) shows results from obtained SuSpect. Gray points in 3(a) and 3(b) are 
consistent with REWSB and LSP neutralino. Green points form a subset of the gray 
and satisfy sparticle mass [35] and B-physics constraints described in Section 4. In 
addition, we require that green points do no worse than the SM in terms of (g — 2) M . 
Brown points form a subset of the green points and satisfy Qh 2 < 1. Panel 3(c) shows 
data collected with Isajet (blue) and SuSpect (red). The red and blue points satisfy 
the constraints imposed on the green points in panels 3(a) and 3(b). 
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5 Higgs Mass Prediction 



In this section we discuss the implicit relationship between the SUSY threshold correc- 
tions to the Yukawa couplings, REWSB and the light CP-even Higgs mass. We begin 
with the SUSY threshold corrections to the top, bottom and tau Yukawa couplings 
which play a crucial role in t-b-r Yukawa coupling unification. In general, the bot- 
tom Yukawa coupling y^ can receive large threshold corrections, while the threshold 
corrections to yt are typically smaller [41]. The scale at which Yukawa coupling unifi- 
cation occurs is identified with Mgut, which is the scale of gauge coupling unification. 
Consider first the unification between ^(M GU t) and ?/ t (M G ut)- The SUSY correc- 
tion to the tau lepton mass is given by 5m T = v cos f38y T . For the large tan/3 values 
of interest here, there is sufficient freedom in the choice of 5y T to achieve yt ~ y T at 
Mqut- This freedom stems from the fact that cos (5 — 1/ tan /3 for large tan /3, and so 
we may choose an appropriate Sy T and tan /3 to give us both the correct r lepton mass 
and y t ~ y T at M GU t- The SUSY contribution to Sy^ has to be carefully monitored 
in order to achieve Yukawa coupling unification ^(Mqut) ~ 2/&(M G ut) ~ Ut(M G u T ). 

We choose the sign of 5yt (i = t, b, r) from the perspective of evolving y± from Mqut 
to Mz- With this choice, 5yb must receive a negative contribution (—0.27 < Syb/yb < 
—0.15) in order to realize Yukawa coupling unification. The leading contribution 
to 5yb arises from the finite gluino and chargino loop corrections, and in our sign 
convention, it is approximately given by [41] 



r finite^ 9z /imgtan/3 y 2 t /iAtan/3 
gb ~ 12vr 2 ml 32vr 2 m? ' 1 ' 

t 

where g% is the strong gauge coupling, rrig is the gluino mass, (m f ~) is the heaviest 
sbottom (stop) mass, and A t is the top trilinear SSB coupling. The logarithmic 
corrections to y^ are positive, which leaves the finite corrections to provide for the 
correct overall negative 5yb in order to realize YU. 

For models with gaugino mass unification or same sign gauginos, the gluino con- 
tribution (first term in Eq.(10)) is positive for /i > 0. Thus, the chargino contribution 
(second term in Eq.(10)) must play an essential role in providing the required negative 
contribution to 5yb- This can be achieved [42, 44] only for 

77ii 6 > Mx/2, m l6 >QTeV and A /m m ~ -2.6. (11) 

One could lower the sparticle mass spectrum by considering opposite sign gaug- 
ino SSB terms [3] which is allowed by the 4-2-2 model [43]. In SO(10) GUT, non- 
universality of SSB gaugino masses with opposite signs can be generated through 
various SO(10) non singlet representations responsible for SUSY breaking [24]. 

In particular, for M 3 < 0, M 2 > and \i > 0, the gluino contribution to 5yb 
has the correct sign to obtain the required 6-quark mass, and furthermore, it is not 
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necessary to have very strong relations among SSB fundamental parameters. Yukawa 
coupling unification in this case is achieved for 

mi 6 > 300 GeV and M 3 (M 2 ) ^ m Wl as well as for M 3 (M 2 ) ^ m ie , (12) 

as opposed to the parameter space given in Eq. (11). This enables us to simulta- 
neously satisfy the requirements of t-b-r YU, neutralino dark matter abundance and 
constraints from (g — 2) M , as well as a variety of other bounds. But for the above 
mentioned cases the relation rn 2 Hd > m 2 Hu was imposed at Mgut- It is well known 
that REWSB cannot be realized if we require m 2 Hd = m 2 Hu at Mgut and, in addition, 
demand exact t-b-r YU (For Quasi YU case [45] one can even have m 2 H = m 2 Hu = m 2 6 
at Mqut). Let's briefly review REWSB condition in light of YU. 
The REWSB minimization condition at tree level requires that 

„ m 2 Hj - m 2 H tan 2 (3 Ml ^ Ml . . 

2 = _JU _Ka > J _ Z 2 Z 13 

F tan 2 (3-1 2 2 v ; 

Here the approximate equality works well for large values for tan (3, which is the case 
for t-b-r YU. Plugging the tree level expression of the CP-odd Higgs mass m 2 A = 
m 2 Hd + rn 2 Hu + 2/i 2 in Eq. (13), we obtain the relation 

m 2 Hd - m 2 Hu >M 2 Z + m\. (14) 

On the other hand, a semi-analytical expression for rn 2 Hd — m 2 Hu at Mz, in terms of 
GUT scale fundamental parameters, has the following form [46] 

m 2 Hd - m 2 Hu w -0.13m 2 6 - 0.26m 2 /2 - 0.04m 1/2 A - 0.01 A 2 , (15) 

which implies that in order to satisfy the condition in Eq. (14), we should require 

M 1/2 > m 16 . (16) 

This requirement clearly contradicts the t-b-r YU condition obtained in Eq. (11) if 
gaugino mass unification is assumed at Mgut- This is the reason why REWSB cannot 
occur if precise Yukawa coupling unification and m 2 H = m 2 H conditions are imposed 
at Mgut- On the other hand, the condition from Eq. (16) can be consistent with 
the condition presented in Eq. (12), with non-universal gaugino masses at Mgut- 
The overlap of conditions from Eqs. (16) and (12) gives very characteristic relations 
among m 16 and gaugino masses. 

We present in Figure 1 the results of the scan over the parameter space listed in 
Eq.(9) in M3 — mie plane. The light blue points are consistent with REWSB and 
Xi LSP. Points in black correspond to 3% or better t-b-r YU (Rtbr < 1.03) for Isajet 
data. For SuSpect data, the black points correspond to 5% or better t-b-r YU. In 
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Section 6 we will discuss the factors that can result in a few percent difference in the 
YU obtained in Isajet and SuSpect. We see that 3% in Isajet's case (5% in SuSpect's 
case ) or better t-b-r YU can occur for M3 slightly heavier than m^, and M3 > 2 
TeV. On the other hand, M3 significantly affects the low scale stop quark mass and 
an approximate semi-analytic expression for for t-b-r YU case is as follows: 

m\ R « 0.27mi 6 + 5.3M 3 2 + 0.4M 2 2 + ... (17) 

Using Eq. (17) we can obtain a rough lower bound for stop quarks, namely mi R > 4 
TeV. 

Next we discuss how the findings presented in Figure 1 can affect the CP-even 
Higgs boson mass calculation. For the actual calculation, both Isajet and SuSpect 
employ a more elaborate calculation procedure. We include the one-loop contribu- 
tions to the CP-even Higgs boson mass [47]: 



[ m h] MSSM ~ M Z COs2 2 ^ ( X 



3 mi \ 3 m\ 



7r 2 v 2 7 An 2 v 2 



t + -X t 
2 * 



M s\ v I A 2 



where 



« = 1741 GeV - ( = 1o H W • Xt = M V ~ nn§) ■ (19) 

Also At = A t — /i cot 0, where A t denotes the stop left and stop right soft mixing 
parameter and Ms = ^jmi L m- tR . Note that one loop radiative corrections to the CP- 
even Higgs mass depend logarithmically on the stop quark mass and linearly on X t . 
These two quantities essentially determine the radiative corrections to the CP-even 
Higgs boson mass. Because of this, it is interesting to present the result from Figure 
1 in terms of X t and Ms- 

In Figure 2 we show the results in the Ms — At /Ms plane. The color coding is 
the same as in Figure 1. From the Ms — At/ Ms plane, we see that the black points 
lie in the interval 5 TeV < Ms < 9 TeV. This interval reduces further if we require 
better YU. This means that for good YU, following Eq.(18) the lightest CP even 
Higgs boson should be relatively heavy owing to the logarithmic dependence on Ms- 

Since the growth or decay of a logarithmic function is slow, the logarithmic de- 
pendence of m/j on Ms nicely explains the shape of colored points in Figure 3(a) and 
3(b). Panel 3(a) shows results obtained from Isajet, while 3(b) shows results obtained 
using Suspect. The gray points in Figs. 3(a) and 3(b) are consistent with REWSB 
and LSP neutralino. The green points form a subset of the gray points and satisfy 
the sparticle mass and B-physics constraints described in Section 4. In addition, we 
require that green points do no worse than the SM in terms of (g — 2) M . The brown 
points form a subset of the green points and satisfy flh 2 < 1. 
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An intriguing feature of Figures 3(a) and 3(b) is that the minima of the distribu- 
tion occurs at the Higgs mass value very close to the observed mass of the SM-like 
Higgs at the LHC. In other words, nearly perfect t-b-r YU prefers the current favored 
value of rrih. We can also understand from Figure 2 why the minima in the Figure 
3(a) and 3(b) have relatively small widths. In Figure 2, from the Ms — A t /M$ plane, 
we see that the ratio A t /Ms lies in the very small interval 0.3< A t /Ms <0.7. On the 
other hand, it is known [48] that the CP even Higgs boson obtains significant contri- 
butions from A t if A t /Ms > 1. We can therefore conclude that there is no significant 
contribution from the finite corrections to the CP even Higgs boson mass if we have 
almost perfect YU, and the Higgs mass is mostly generated from the logarithmic 
corrections. This is why the minima in Figures 3(a) and 3(b) are not wide. 

In Figures 3(c) we show the overlap of the data from Isajet 7.84 and SuSpect 
2.41. All the points shown are consistent with REWSB, LSP neutralino, and satisfy 
the sparticle mass and B-physics constraints described in Section 4. The blue (red) 
points show results obtained from Isajet (SuSpect). We can see that the minima 
of the data distribution from Isajet predicts a Higgs mass rrih ~ 124 GeV, whereas 
SuSpect predicts m h ~ 126 GeV. For lower values of Rtb T , the results from the two 
packages overlap around rrih ~ 125 GeV. This observation makes the predicted value 
of Higgs mass close to 125 GeV obtained from t-b-r YU condition quite reliable. 

6 Higgs and Sparticle Spectroscopy From Isajet 
and SuSpect 

In this section we compare the allowed parameter spaces and the sparticle spec- 
troscopy obtained from Isajet and SuSpect. A comparison of the uncertainties in the 
sparticle spectroscopy of different packages was done in [32], and, in particular, the 
threshold effects were compared in [33]. The approach used by the two programs 
is very similar in that both use two loop RGE running for the gauge and Yukawa 
couplings. But there are some factors that can lead to numerical differences of a 
few percent. These include the scale at which the sparticles are integrated out of 
the theory. In Isajet the SSB parameters are extracted from RGE running at their 
respective mass scales m; = mj(mj), whereas in SuSpect these parameters are ex- 
tracted at Mewsb(= ■J m t-J n h)- SuSpect uses a s in the DR scheme, whereas Isajet 
uses the MS value. Another source of discrepancy can be the use of bottom pole 
mass by SuSpect rrib = rrib(Mz), whereas Isajet uses the mass at the SUSY scale 
m b = m h (M S usY)- The default guess of the sparticle masses at the beginning of the 
RG evolution process is different in Isajet and SuSpect. SuSpect assumes that the 
full set of MSSM RGEs are valid between M z and M GUT . This is also true for Isajet, 
except that it employs one-loop step-beta functions for gauge and Yukawa couplings. 
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Figure 4: Plots in the M1/2 — ttiiq, Aq/ttiiq — m 16 , m 10 — m 16 and [L — m 16 planes. All 
the points shown are consistent with REWSB, LSP neutralino and satisfy sparticle 
mass [35] and B- physics constraints described in Section 4. We also require that 
the points do no worse than the SM in terms of (g — 2) M . In addition the satisfy the 
Higgs mass range 122 GeV < rrih < 128 GeV and R t b T < 1.2. The purple points show 
results obtained from SuSpect and yellow points is the data collected using Isajet. 
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In Figure 4 we show our results in the My 2 — mm, ^o/ m w — m w, m w ~ m w 
and \i — rriiQ planes. All of the points shown are consistent with REWSB, LSP 
neutralino and satisfy sparticle mass and B-physics constraints described in Section 
4. In addition, the points shown satisfy the condition 122 GeV < < 128 GeV 
and R t br < 1-2. The purple points show results obtained from SuSpect, whereas the 
yellow points correspond to the data collected using Isajet. 

We observe that there are some small but notable differences in the results ob- 
tained from the two packages. It is interesting to note that for the regions where 
the points are more dense there is good overlap between the two packages. In Fig- 
ure 4 under the very dense yellow points there also lie dense purple points. From 
Mi/2 — miQ, A /miQ — miQ, m 10 — m 16 panels, we see that the best agreement between 
solutions obtained from Isajet and SuSpect occurs for m 16 < 6 TeV, M 1 / 2 < 3 TeV 
and m w < 4 TeV. 

In the fi — rriiQ plane of Figure 4 we see that the solutions from SuSpect with 20% 
or better t-b-r YU have lower values of fi compared to the /i values from Isajet. To 
exemplify this further, in Figure 5 we show a plot in R t b T — fi plane. We see that 
requiring 20% or better t-b-r YU allows solutions from SuSpect with \i <m 200 GeV, 
while requiring the same for Isajet yields \i > 1 TeV. On the other hand, requiring 
YU better then 5% leads to similar limits on the values of [i from Isajet and SuSpect, 
i.e., n > 2 TeV. It was noted in [34] that the two codes can differ notably in regions 
with low fi and tan /3 values. This difference can stem from the factors previously 
mentioned, and which may have important implications for natural SUSY. 

We observe from Figure 5 that SuSpect does not yield YU better than 5%, with 
a minimum value of Rtbr ~ 1.05. Isajet, however, predicts even better YU with the 
minimum value of Rtbr ~ 1.02. This few percent difference could be due to the way 
threshold effects are evaluated by the two codes [33]. The R t b T — m 16 plane shows 
that in order to have YU better than 5% the results from Isajet require 600 GeV 
< m 16 < 2.5 TeV. 

In Figure 6 we show results in the R t b T — mA, Rtbr — m^, Rtbr — Rtbr — m§, 
Rtbr — rn-x anc ^ ^ tbr ~ m x ± pl anes - The color coding is the same as in Figure 3. 
The data collected with Isajet is used to make the plots in this figure. All panels in 
Figure 6 indicate that the model predicts relatively narrow ranges for the sparticle 
masses corresponding to the best t-b-r YU. The sparticles are heavy enough to evade 
observation at current LHC energies, but a signal would hopefully be observed during 
the 14 TeV LHC run. 

From the Rtb T — m>A plane we see that just from the REWSB condition, t-b-r YU 
better than 5% predicts that we cannot have a very light CP-odd Higgs boson (A). 
Similar bounds apply for the heavier CP even H boson and charged H bosons, since, 
in the so-called decoupling limit, m\ ^> Mf , we have m 2 H ~ m\ and m#± ~ m\+M^ v 
[49] , where Mw stands for the PU-boson mass. The leading decay modes for the heavy 



15 




Figure 5: Plots in the R t b T — A*, Rtbr — m w planes. Color coding is the same as in 
Figure 3(c). 



Higgs bosons are H, A — >■ bb and H,A^r t + t~ . The heavy Higgs production cross 
section and branching ratios depends on tan (3 and its mass. The MSSM sparticle 
dependence appears through gluino-squark, Higgsino-squark, bino-sfermion and wino- 
sfermion loops. The bound on tua in our scenario is more relaxed since we have non 
universal gauginos with opposite signs for M 2 and M 3 at the GUT scale. We also 
find that good t-b-r YU requires M2/M3 > 1.5. This will alter the gluino-squark and 
wino-sfermion loop contributions to the heavy Higgs production cross section and 
branching ratios [50] . 

Applying all the collider and B-physics constraints we obtain a lower bound for 
rriA which is very close to the value corresponding to best YU. Restricting to 5% 
or better unification and including the constraints presented in Section 4, we obtain 
the bound 400 GeV< Ma < 1 TeV. The lower bound is very close to the current 
experimental limit [10] obtained from the GUT scale gaugino unification condition 
and can be further tested in near future. 

In the R t br — to^ plane we can observe that the preferred values for the stau 
lepton mass from the point of view of good YU (Rtbr < 1-05) is in the interval 500 
GeV < m? 1 < 1.5 TeV. The search for a stau in this mass range is challenging at 
the LHC. In this model, the stau is the NLSP, and this can yield the correct relic 
abundance through neutralino-stau coannihilation. 

The lightest colored sparticle in this scenario is one of the stops. The preferred 
mass as seen from the R t b T — Tn^ plane, is around 3-4 TeV. For the gluino, the mass 
according to the R t b T — plane is around 5-6 TeV. In principle they can be found at 
the LHC. We can also see from Figure 6 that the lightest neutralino, for R t b T < 1.05, 
is around 500 GeV, and the preferred value from the point of view YU is ~ 700 GeV. 
The model also predicts the charginos to be heavier than 2 TeV. 

Since the pseudoscalar A boson is being searched for at the LHC, we present in 
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Figure 6: Plot in the R thr - m A} R tbr - m fl , R tbT - m iv R tbr - m tadeg , R t b T - m^o 
and Rtbr — Tn%± planes. The color coding is the same as in Figure 3. 
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Figure 7: Plot in the R tbr — m A planes. Color coding is the same as in Figure 3(c). 

Figure 7 the combined results from Isajet and SuSpect. The color coding is the same 
as in Figure 6, and we see that the agreement between the two programs is quite 
satisfying. 

In Figure 8 we show the implication of our analysis for direct detection of dark 
matter. Plots are shown in the ogi — m x\ an< ^ °sd — m x ^ planes. The gray points in 
the figure are consistent with REWSB and LSP neutralino. The green points form 
a subset of the gray points and satisfy the constraints described in Section 4. The 
brown points form a subset of the green points and satisfy Qh 2 < 1, and the orange 
points form a subset of the brown points with 10% or better YU (Rtbr < 1-1). The left 
panel shows the current and future bounds from CDMS as black (solid and dashed) 
lines, and as red (solid and dotted) lines for the Xenon experiment. The right panel 
also shows the current bounds from Super K (solid red line) and IceCube (solid black 
line), and future reach of IceCube DeepCore (dotted black line). 

We can see that the parameter space of this model representing neutralino-stau 
coannihilation can be tested with these experiments. However models with YU better 
than 10% yield tiny cross sections which are well below the sensitivity of these exper- 
iments. We also observe that good YU predicts a heavy neutralino (tu^o > 400 GeV). 
It is, however, interesting that these orange points also predict a 125 GeV Higgs mass 
as seen in Figures 3(a) and 3(b). Therefore, t-b-r YU not only predicts a 125 GeV 
Higgs but also a relatively heavy dark matter LSP which coannhilates with the stau 
to yield the correct relic abundance, and also yields tiny cross sections well below the 
sensitivity of current experiments. 

In Table 1 we present four benchmark points with good YU and Higgs mass ~ 125 
GeV. The points shown also satisfy the constraints described in Section 4. Points 1 
and 2 represent solutions that yield the best YU in Isajet and SuSpect. As described 
earlier, Isajet yields YU as good as ~ 2%, whereas in SuSpect it is ~ 5%. Point 3 
depicts stau coannihilation in addition to a 124 GeV Higgs and R t b T = 1-03. Point 4 
shows that good YU can be attained with the sfermions and Higgs nearly degenerate 
at M GUT , i.e., m 16 ^ m 1Q . 
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Figure 8: Plots in the o^i — m^o and o"sd — m x ° planes. The cross sections are 
calculated using Isajet. Points shown in gray are consistent with REWSB and LSP 
neutralino. Green points form a subset of the gray and satisfy sparticle mass [35] 
and B-physics constraints described in Section 4. In addition, we require that green 
points do no worse than the SM in terms of (g — 2) M . Brown points form a subset of 
the green points and satisfy Qh 2 < 1. The orange points are a subset of the brown 
points and satisfy R t b T < 11 In the crgi - m x o plane, the current and future bounds 
from the CDMS experiment are represented as black (solid and dashed) lines and as 
red (solid and dotted) lines for the Xenon experiment. The right panel shows the 
c"sd - m x\ plane with the current bounds from Super K (solid red line) and IceCube 
(solid black line) and future reach of IceCube DeepCore (dotted black line). 




Figure 9: Plots in the R t b T -^EW and R t b T -^HS planes. Color coding is the same as 
in Figure 3. 
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7 Fine tuning constraints for little hierarchy 



The latest (7.84) version of ISAJET [29] calculates the fine-tuning conditions related 
to the little hierarchy problem at Electro Weak (EW) scale and at the GUT scale 
(HS). We will briefly describe these parameters in this section. 

After including the one-loop effective potential contributions to the tree level 
MSSM Higgs potential, the Z boson mass is given by the following relation: 



Mf _ {m\ d + Sj) - {m\ u + Eg tan 2 g 
2 " tan 2 (3 - 1 



A* • (20) 



The E's stand for the contributions coming from the one-loop effective potential (For 
more details see ref. [51]). All parameters in Eq. (20) are defined at EW scale. 

In order to measure the EW scale fine-tuning condition associated with the little 
hierarchy problem, the following definitions are used [51]: 

C^ = |m^/(tan 2 /3-l)|, C Hu = | - m 2 Hu tan 2 /3/(tan 2 (3 - 1)|, C M = |-/i 2 |, (21) 

with each C 2 u,d n less than some characteristic value of order Mf . Here, i labels the 

u,d y t 

SM and supersymmetric particles that contribute to the one-loop Higgs potential. 
For the fine-tuning condition we have 

A EW = max(C l )/(Ml/2). (22) 

Note that Eq. (22) defines the fine-tuning condition at EW scale without addressing 
the question of the origin of the parameters that are involved. 

In most SUSY breaking scenarios the parameters in Eq. (20) are defined at a 
scale higher than Mew- in order to fully address the fine-tuning condition we need 
to check the relations among the parameters involved in Eq. (20) at high scale. We 
relate the parameters at low and high scales as follows: 

m k >d = m kA M Hs) + Srn 2 Had , ^ = fi 2 (M HS ) + 5fi 2 . (23) 

Here m 2 Hu d (Mns) and /i 2 (M#s) are the corresponding parameters renormalized at 
the high scale, and 5m 2 Hud , 5{i 2 measure how the given parameter is changed due to 
renormalization group evolution (RGE). Eq. (20) can be re-expressed in the form 

m| = (m 2 Hd (M HS ) + 5m 2 Hd + Sj) - (m 2 Hu (M HS ) + 5m 2 Hu + Eg tan 2 g 
2 tan 2 /3 - 1 

- ( l x 2 (M HS ) + 5 l i 2 ) . (24) 

Following ref. [51], we introduce the parameters: 

B Hd = \m 2 Hd (M HS )/(t^ 2 f3 - l)\,B SHd ee \5m 2 H J(t^ 2 f3 - 1)|, 
B Hu ee | - m 2 Hu (M HS ) tan 2 (3/ (tan 2 (3 - 1)\,B^ = \fi 2 (M HS )\, 

B 5Hu ee | - 5m 2 Hu tan 2 /3/(tan 2 /3 - 1)|,£^ = (25) 
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and the high scale fine-tuning measure Ahs is defined to be 

A HS = max(Si)/(M|/2). 



(26) 



The current experimental bound on the chargino mass {m^ > 103 GeV) [35] 
indicates that either A EW or A HS cannot be less than 1. The quantities A EW and 
Ahs measure the sensitivity of the Z-boson mass to the parameters defined in Eqs. 
(21) and (25), such that (100/A W )% ((100/A HS )%) is the degree of fine-tuning at 
the corresponding scale. 

Based on the definition of high and low scale fine tuning described above we show 
results in the R t b T — A EW and R t b T — A HS planes (using Isajet) in Figure 9. We see 
that the low scale little hierarchy problem becomes more severe for the t-b-r YU case 
compared to what we have in the constrained MSSM (CMSSM), but high scale fine 
tuning is at the same level as in the CMSSM [52]. For t-b-r YU of around 5%, the 
EW fine tuning parameter A E w ~ 800 and the HS fine tuning parameter is also 
Ahs ~ 800. As mentioned above, the fine tuning condition has to be scale invariant 
which means that cancellation between parameters at a particular scale cannot be 
more severe compared to same conditions at another scale. Based on this assumption 
the little hierarchy problem in this model remains the same as we have when gaugino 
universality is assumed in the theory [51, 52]. 

8 Conclusion 

We have demonstrated how t-b-r YU is consistent with a 125 GeV Higgs. Our analysis 
is an extension of the analysis in ref. [6] in several ways. We have highlighted the 
effects of threshold corrections on the bottom quark Yukawa coupling in this model, 
and discussed the implicit relationship between these corrections and the Higgs mass. 
We showed that for YU better than ~ 5%, M 3 > m 16 at M GUT , with M 3 > 2 TeV. 
This, in turn, leads to a heavy stop quark, m^ R > 4 TeV. The dominant contribution 
to the Higgs mass arises from the logarithmic dependence of rrih on the stop quark 
mass. This leads to the prediction rrih ~ 125 GeV, consistent with t-b-r YU better 
than 5%. 

We also compared our results from two different packages, namely Isajet and 
SuSpect. We found good agreement between the two codes with only a few percent 
difference between the calculations. One important difference is that Isajet allows YU 
better than 2%, whereas SuSpect has, at best, 5% YU. Another notable difference 
is that SuSpect allows for much smaller values of the Higgs mixing parameter \i. 
This can have implications for natural SUSY since smaller ii values are preferred in 
resolving the little hierarchy problem. 

The two codes also agree well in their predictions of sparticle masses. We find 
that insisting on YU better than 5% implies that the sparticles are heavy enough 
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Tcq ipf 

lOd 1 v3 u 


'"10 

mi 6 
Mi 
M 2 

/I n ITT) 1 c 
-<M)/ '"'16 

tan /3 
mt 


4.19 x 10 2 
2.13 x 10 3 
1.89 x 10 3 
5.67 x 10 3 
—3 78 x 10 3 
2.39 
47.18 
174.2 


3 82 x 1 3 
2.69 x 10 3 
2.00 x 10 3 
6.00 x 10 3 
—4 00 x 10 3 
1.37 
48.05 
173.1 


4.49 x 10 2 
1.91 x 10 3 
1.78 x 10 3 
5.35 x 10 3 
—3 57 x 10 3 
0.03 
47.93 
174.2 


1.94 x 10 3 
2.00 x 10 3 
1.51 x 10 3 
4.53 x 10 3 
—3 02 x 10 3 
1.56 
47.46 
173.1 


H 


3729 


1935 


2913 


2526 


m h 

m H 
m A 
m H ± 


125 
747 
742 
753 


126 
491 
491 
500 


124 
572 
568 
580 


123 
558 
554 
567 


X-1,2 


895, 3739 


955, 1935 


848, 2932 


709, 2540 


X-3,4. 
*1,2 

m g 


3742, 4822 


1936, 5043 


2935,4562 


2543, 3849 


3789, 4774 


1934, 5043 


2978,4516 


2579, 3809 


7694 


7673 


7266 


6239 




7667, 6824 
5331,6560 


8112,7245 
5604, 6839 


7219,6415 
5239,6367 


6295, 5635 
4390, 5370 


m d L , R 
mi 

01,2 


7668,6814 
5553, 6526 


8112,7236 
5870, 6870 


7220, 6406 
5434, 6333 


6296, 5628 
4591,5341 




4148 
3898 


4590 
4234 


3870 
3641 


3487 
3243 


m e L , R 

m fl2 


4153,2234 
1094, 3875 


4590, 2780 
1140,4235 


3875, 2009 
881,3620 


3491, 2068 
1061,3225 


A(<? - 2)„ 
o"5/(pb) 
vsd (pb) 

&CDMh 2 


3.11 x 10- 11 
1.59 x 10- 11 
4.69 x 10~ 10 
6.5 


3.36 x 10- 11 
1.29 x 10~ 9 
1.35 x 10- 9 
2.8 


3.71 x 10- 11 
7.08 x 10- 11 
1.60 x 10~ 9 
0.8 


4.97 x 10~ n 
1.00 x 10~ 10 
2.89 x 10" 9 
4.0 


Rtbr 


1.02 


1.05 


1.03 


1.04 



Table 1: Benchmark points with good Yukawa unification. All the masses are in 
units of GeV. Point 1, 3 and 4 are generated using Isajet 7.84 whereas point 2 is 
from SuSpect 2.41. Point 1 and 2 demonstrates how a small value of Rtb T yields a 
Higgs mass ~ 125 GeV. Point 3 exhibits stau coannihilation and has a small R t b T 
that agrees with Qh 2 < 1. Point 4 has m 16 ~ m 10 with good YU. 
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to evade observation at the current LHC energies, but may be observed during the 
14 TeV LHC run. Furthermore, we showed that t-b-r YU predicts a light CP-odd 
Higgs boson (A). Restricting to 5% and better YU yields the following bound on the 
pseudoscalar mass, 400 GeV< Ma < 1 TeV. Similarly, the bounds on other sparticle 
masses are rrig > 4 TeV, m f > 500 GeV and m^± > 2 TeV for YU 5% or better. 

Finally, we also tested the implications of YU for direct detection of dark matter. 
We found that stau-neutralino coannihilation can lead to the correct dark matter 
relic abundance. Moreover, insisting on YU better than 10% implies a heavy dark 
matter candidate (m^o > 400 GeV). The neutralino-nucleon cross sections are found 
to be well below the current sensitivity of direct detection experiments. 
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